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1 Introduction 



Physics of QCD jets is one of the most important testing grounds for the theory of strong inter- 
actions. The perturbative evolution of a quark-gluon jet is at present well understood (see, e.g., 
The known details include in particular the famous colour coherence phenomenon which leads 
to a probabilistic picture of timelike jet evolution where the daughter partons are emitted into a 
gradually shrinking cones. The nonperturbative aspects of a jet evolution are much less understood. 
From the experience of a sum rule approach applied to the analysis of heavy resonance properties 
|^] it is clear, that when the virtuality of a daughter parton is of order of 1 GeV one has to take 
into account the nonperturbative corrections due to the nonperturbative vacuum fields giving rise 
to the QCD vacuum condensates. Another situation in which one deals with dissipative effects is 
a jet propagation in, e.g., nuclear medium. Although physically these situations are very different, 
we expect the corresponding theoretical formalism to be the same. 

In both cases the new interaction vertices lead to the appearance of new dimensionful param- 
eters. This means in turn that a scaling description of jet evolution where the interaction vertices 
are determined only by dimensionless quantities (energy or virtuality ratios) is no longer valid. 
Physically what happens is a beginning of a string formation taking energy from the perturbative 
component and converting it to nonperturbative degrees of freedom. The first model description 
of a jet evolution taking into account the nonperturbative energy loss was proposed long ago by 
Dremin Q. It was based on the analogy with the physics oh high energy electromagnetic showers 
in the medium, where apart from a scale-invariant evolution due to a photon bremmstrahlung and 
pair creation there appears a scale-noninvariant energy dissipation due to the ionization of the 
atoms in the medium j|. Later the corresponding modified evolution equations were analytically 
solved in ||, where the expressions for the parton multiplicities in quark and gluon jets and the 
energy loss by the perturbative component were computed. A discussion of this approach and also 
of the latest Monte-Carlo calculations, || where the nonperturbative component was described by 
QCD effective lagrangeans can be found in the recent review jjj. A serious drawback of Dremin 
equations || is, however, a disregard of the kinematical and color interference effects in the jet 
evolution which are known to be of crucial importance for describing the jet characteristics. The 
main goal of the present paper is to introduce a formalism providing a possibility of an analyti- 
cal treatment of the jet evolution taking into account both color interference and nonperturbative 
energy loss. 

Before writing down the modified evolution equations let us recall, that color coherence in QCD 
jets can be described using various approximations. The simplest is a leading logarithmic one 
(DLA, see [|]]). The calculations in this approximation are relatively simple, but the energy of the 
jet is not conserved. This is due to the eikonal description used in this approximation, when the 
energy loss of the projectile is considered to be negligible. A more refined analysis is made within 
a modified leading logarithmic approximation (MLLA), where the total energy is conserved. The 
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Dremin equations H correspond to a case, where the energy is perturbatively conserved. Thus to 
calculate an absolute value of the dissipative energy loss one has to modify an MLLA formalism. 
The corresponding equations can be written (see Appendix A), but unfortunately we were not able 
to find their solutions. Thus in what follows we shall concentrate on a simpler DLA case, where 
the energy nonconservation is built in already at the level of perturbation theory. In this case one 
can not calculate the absolute energy loss, but the relative one is fairly well-defined. 



2 DLA Equation for the Generating Functional 

Let us turn to a description of a formalism which generalizes the known effective methods of 
describing the QCD jet evolution including color coherence effect [jjj by providing a possibility of 
considering the new possible nonperturbative sources of energy dissipation. 

Technically it is most convenient to use a generating functional. G({u}) from which one can 
calculate both an exclusive cross section 

<=(nft4) c,,,|| L>' (1) 



and an inclusive one 



N 



\i=l 



In the following we shall use the notations dK = 7ox 2 fc^ ^ or ^ ne usual DLA integration measure 
and T(p, 9) for the possible gluon emission phase space domain, where p is an initial gluon energy 
and 9 is an opening angle of a jet. Let us remind, that in DLA approximation the generating 
functional G(p, 9; {u}) satisfies the following master equation 



G(p,9;{u}) =u(p)e- w( -P^ + 

- wM+wip ' ek) G(p,9 k ;{u})G(k,9 k ;{u}), (3) 



dk d 2 k±_ 2Cya s 
r( P ,0) k ki 2 vr 



where 



w(p,9) = / dK 
Jr(p,6) 



is a total probability of an emission of a gluon having an energy p into the cone having an opening 
angle 9, and exp(— w(p, 9)) is a formfactor giving a probability of a non-emission of the same gluon. 
Then a sandard procedure leads to the following equation on the generating functional: 

G(p,9;{u})=u(p)exp( f dK [G(k, 9 k ;{«})- l]j . (4) 

Let us now discuss the possible nonperturbative modifications of the equations for the generating 
functional. The aim is to describe a partial convertion of the perturbative degrees of freedom 



3 



(gluons) and their energy into the nonperturbative ones by the process presumably analogous to 
the ionization in the case of electromagnetic showers discussed in the Introduction. Physically 
the picture for the jet is that the dissipation pumps some energy from the cone corresponding to 
the perturbative jet evolution reducing the gluon multiplicity and energy inside it. Thus we have 
some additional nonperturbative vertex coupling the perturbative gluons with the nonperturbative 
vacuum or nuclear medium. The presence of such a coupling can be related to reducing a probability 
of an emission of a perturbative gluon by taking away a certain portion of its energy before the 
perturbative emission. Techically it is convenient to introduce a following modification of the gluon 
emission probability: 

w(p,6-0) = [ dK, (5) 

where the constant (3 describes the additional damping of a perturbative gluon emission due to a 
dissipative interaction with the nonperturbative vacuum fluctuations or nuclear medium introduced 
by compressing a possible phase space domain for the gluon emission from Y{p, 9) to T(p — (3,9). 
The corresponding modification of the master equation (^) reads 

G(p,9;{u}) =u(p)e- w (P' 9 ^ + 

dk tfkL 2( V ° s - W <P*fO+°<P**> G(pt 9k ; {u})G{k, 9 k ; {«}), (6) 



h( P -f3,e) k k± 2 it 
Then instead of the Eq. (Q) we get 

G{p, 9; {u}) = u(p) exp ( [ dK [G(k,9 k ; {u}) - 1] ] , (7) 

where the initial conditions for this equation read G(p,9; {u})\ u =i = 1. 

We suggest that /3<j) for all physically meaningful values of energy and in the following we 
shall consider the terms of the first order in the small parameter e = (3 /p. In Appendix B we 
illustrate the thus arising perturbation theory by deriving the differential equation on the leading 
correction to the generating functional. 

3 Mean multiplicity 

Let us begin this section by looking at the simplest jet characteristic, a mean multiplicity of partons 
in a jet (in the following we are considering only gluons). Following the standard steps we obtain 
a following integral equation on the mean multiplicity n(p, 9): 

n(p,9) = 1+1 dKn{k,9 k ). (8) 

In the following it will be convenient to rewrite the integration over the phase space as 

/ dK=f f fpr7o 2 =/^ fdy'il (9) 

J T{p- 0,6) JT{p-l3fi) k lTXk\ J J 
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where 

. pe , ko k ko 

y = m—,y =ln— = In— . 
Vo Wo 

Working in the first order in e we can write 

n(p,6) = n (y) +efi 1 {y). 

The equations for the functions no(y) and ni(y) take the form 

y 

no(y) = l + j dy'^n (y') (10) 
o 

and 

y y 
ni(v) = J dy'^fniy') (e^' - 1) - | c^To 2 ^') (H) 
o o 
The first one has a well known solution 

no(y) = ch(7 y), 
and the second one can be rewritten in the differential form 

n'iiy) - n[(y) = 7o™i(y) + 7o( ch (7oy) - 7osh(7 y)) (12) 
with the initial conditions n(0) = 0,n'(0) = — 7q- The solution reads 

\ 2 e Xiy - Xie X2V 

ni = 7o cn (7oy) - 7osh(7o?/) + +70 , — • (13) 

V 1 + 4 7o 2 

In the limit of (y 1, 70 <C 1) we have 



BM ..™ (1 _ rB _^ ) = .^j ^- 70 ^- 2 ^j (14) 

The above result is in accord with the expectation based on physical reasoning. Namely, the 
dissipation should lead to a decrease in perturbative multiplicity. Let us note, that from Eq. ( |l4|) 
it is clear, that in order for perturbation theory to be applicable one should have (3 <C p/70 and 
(3 <C Qo/lo- This means that apart from the expected expansion parameter e there appears a new 
one e' = An interesting feature of the above answer is a multiplicative combination of the 

dissipation scale (3 and a perturbative coupling constant hidden in 70. 
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4 Energy distrubution of particles in a jet 



Let us now turn to a computation of a one-particle energy distribution function D(k, 9). It can be 
obtained from the generating functional in the following way: 

D(k, 9) = k^^G(p, 6; W)[ =i , (15) 

Making use of (0) and (|) we get 

l-e y 

D(l,y) = 8(1) + / dl' [ dy'^Dil^y'), (16) 





where I = ln(p/k) = ln(l/x) and V = ln(k' /k). The integration is performed over k' and 9y. Let us 
stress that in the above equation (|i~6| ) the distrubution function D is also (3 - dependent. Expanding 
it in the small parameter e 

D(l,y) = D (l,y)+eD 1 (l,y) 

we rewrite ( |l6| ) in the form 

l-e y 

D (l,y)+eD 1 (l,y) = 8(l)+ J dl' J ' dy' 7 i(D (l,y) + ee l ~ v ' D x {l,y)). (17) 

o o 

Considering the zero and first order in e we have 

I y 

D (l,y) = 8(1) + j dl' J dy'~f 2 D (l',y'). (18) 


i y 

3,(1, y) = J dl' J dy'^e 1 - 1 ' 'D 1 (l',y') - J" ' dy'^D^y'). (19) 
o o 

The solution of equation ( |18D reads 

D (l,y) = 8(1) +^^h(2 l0 Vyl) (20) 
and that for D\ is (for details see Appendix C): 

Di(l,y) = -Vll8(l) - il(e l + l)|l 2 (2 7 o^/) + 7o (f ^ & ~ 1)J 3 (2 7 oa/^)- (21) 

The plot of the resulting distribution for the jet energy 20 GeV is given at Fig. 1. We see that, 
as anticipated, the effect of dissipation shows itself through the noticeable reduction of distribu- 
tion function. Let us also mention that here we also have a multiplicative dependence on both 
perturbative and nonperturbative factors. 
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5 Conclusions 

In this paper we have proposed a phenomenological nonperturbative modification of the equations 
describing the evolution of QCD jets and taking into account the crucial feature of the colour coher- 
ence of the QCD cascades by accounting for nonperturbative energy dissipation. The corresponding 
modification of the MLLA and DLA formalism exploiting the generating functional was proposed. 
The calculation of simplest jet characteristics such as mean multiplicity and energy spectrum of 
particles in a jet in a DLA approximation has demonstrated an expected decrease in multiplicity 
and corresponding changes in the energy distribution. 

An interesting feature of the result is an unusual perturbative damping of the introduced non- 
perturbative energy dissipation appearing in the multiplicative dependance on some power of the 
QCD coupling constant times the dissipative scale. It is tempting to relate this feature with the 
successes of the soft blanshing hypothesis, where it is assumed that the nonperturbative effects are 
not crucially essential for the jet characteristics. It is interesting to see, whether such damping is 
present within a more realistic MLLA description. Work in this direction is in progress. 
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Appendix A 



Let us start with the proposed general modification of the equation for the generating functional 
(we are closely following the standard notations of Q): 



G A (p,0;{u}) = u(p) e~ WAM + dK 7 g(jfc ± ) e -^W)+»A(pA) 

Jr( P ,e) 

\k1 c £) G B (k,9 k ;{u}) Gc(p-k,9 k ;{u})- 



P 



(A.l) 



Lf(k,p) G B (k,9 k ;{u}) 



Here L^(k,p) is a phenomenologically introduced vertex responsible for the nonperturbative energy 
dissipation. It is convenient to rewrite the above equation in the form: 



G A {p,0]{u}) e WAM = u(p)+ dK-y$(k x ) e «*(p,e*) * 

Jr( P ,8) 

.2 " p 

-L%(k,p) G B (k,6 k ;{u}) 



K% C (-) G B (k,6 k ;{u}) G c (p-k,0 k ;{u})- 



(A.2) 



This equation has to be supplied by the boundary condition Ga = 1 at u = 1. This gives an 
equation of the formfactor wa(p,0)'- 



1 + 



r(p,0) 



dK^(k ± ) e WA{p ' 6k) 



n n P r 



where z = |. Standard procedure leads us to the following master equation: 



a in e 



dz 7q (fcj 



\K% c (z) {G B (zp,B) Gc((l-z)p,e)-J2G A (p,e) 



B,C 



-L B A (z.p,p)\G B (zp,6) -^G A (P,9) 



(A.3) 



Let us now consider a pure qluodynamics and a specific form of the nonperturbative vertex 
&_ 

p 



L^(z ■ P,p) = v 6(1 — ^ — z). Then the master equation takes the form 



dG 



i 



dz 7o( fc ±) 



K{z) {G(zp, 9) G((l - z)p, 6) - G(p, 6)} 



2 

v^(k±){G(p-(3',e) -G(p,0) } 
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In the interesting case of (3 <C p we get 

G( P -f3',6 



G(jp,0) 



-ft' 



and 



,dG 
dp 



dG 
d\n9 



-K(z) {G(zp, 6) G((l - z)p, 9) - G(p, 6)} 



+ «f. (A.4) 



The above equation contains a characteristic derivative over energy (the last term). This term 
is responsible for the energy dissipation and has previously been used for describing the dissipative 
energy losses in the electromagnetic showers B and QCD cascades M, M. 



Appendix B 

In this Appendix we derive an equation for the first correction to the generating functional. Let us 
start with the basic equation 



G(p, 9; {u}) = u(p) exp 



(f dK[G(k,d k ;{u})-l}) 



(B.l) 



with the initial condition G(p, 9; {u})\ u= i = 1. Expanding the generating functional we get in the 
zeorth and first orders in s: 



G(p, 9; {u}) = G (y; u) + eG^y; u), 

where y = In(p9/Q ), G (y; 1) = l,Gi(y; 1) = and e = (3/p. 
An argument of the exponent in flB.l| ) is 

v y 

J dy' 1 2 (y-y')(G (y';u)-l)-e J dy'^(G (y';u - 1))- 
o o 



(B.2) 



-e j dy'rtiey-y'^iy'-u). 
o 

After expanding this exponent to the first order in e we arrive at 



G (y;u) = uexp 



dy'ioiv ~ y')(Go{y';u) - 1) 



(B.3) 



(B.4) 
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and 

Gi{y-u) = G (y;u) 
or in a simpler form 



dy'll{e»-y' - l)Gi(i/;u) - / dy' ^(G (y';u) - 1) 



and 



G 



G '(y) = G (y) J d y '^(G (y') - 1) 
o 

y 

(y) = G (y) J dy'^-y' - 1)^(2/') - G '(y) 



The further simplification can be achieved by introducing a new function H (y) by 

y 

H(y) = J dy'( e y-y' - l)Gi(i/) Gi'(y) = £T"(y) - H' (y) 



(B.5) 



(B.6) 



(B.7) 



(Bi 



After the simple calculation we obtain the following differential equation on the leading correction 
to the generating functional 



H"(y) - H'(y) = ^G (y)H(y) - G '(y) 



(B.9) 



The solutions of the above equation could be used as a starting point for the analysis of the 
dissipative corrections to the scaling behaviour of the moments and, in general, to the KNO property 
of the leading approximation. 



Appendix C 

In this Appendix we shall give some details on the computation of the correction to the distribution 
function D(l,y). We start with the equaitons Eqs. (18), (19): 

D(l,y) = D (l,y) + eD 1 (l,y) 
( y 

'.,2 fS til „/\ 



D (l, y) = 8(1) +J dl'J dy'^D (l', y'). (C.l) 
o o 

l y 

£i(Z,y) = J dl'J dy'^e 1 - 1 ' Si(Z',y') - jT dy' 7 2 D (l,y'). (C.2) 



o o 



The solution of equation (18) is performed by substituting 

D (l,y) = 5(l)+J2C^ n l m y n , 
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giving 



and 



m,n — °m,m+l 



2(m+l) 

m!(m + 1)! 



A)(*,l/) = *(0+W yA(27o VI/O 
Let us look for the solution for Z?i in the form 

The solution for reads 



C (i) 

^ rn r) 



2n n-2 
m+n To /om— n+2 



V L ) n\m\ ^ u m-p 
p=0 



7o <m+1) 
+■ <^m+2,n I — ° m \ n \ 







Thus we obtain the correction for the energy distribution 
Di(l,y) = -yil5{l) - e l 



n < m + 2 
otherwise 



2n "-2 

H2 10 Vyi) + E(-) m+n ^ri E <Cr£ +2 ' T v 

* — ' m\n\ * — ' F 



m,n 



p=0 



which can finally be be written in a compact analytical form Eq. (|21|): 

Si(*,lO = -Vllm ~7oV + l)f / 2 (2 7 o^) + 7o 3 (f ) 3/2 (e* - l)/ 3 (2 7 oV^)- 



11 



